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Abstract 

The next-to-leading order analysis of the cross section of hadronic gluino pair 
production close to threshold is presented. Within the framework of non-relativistic 
QCD a significant enhancement compared to fixed order perturbation theory is 
observed which originates from the characteristic remnant of the IS peak below 
the nominal pair threshold. This enhancement is similar to the corresponding one 
for top production. However, as a consequence of the larger colour factor of the 
QCD potential the effect is significantly enhanced. The analysis includes all colour 
configurations of S-wave gluino pairs, i.e. singlet, symmetric and antisymmetric 
octet, decuplet and twenty-seven representation. Matching coefficients involving real 
and virtual radiation are separately evaluated for all colour and spin configurations 
and initial states. We concentrate on the case of gluino decay rates comparable 
to the gluino binding energy. The non-relativistic dynamics of the gluino pair is 
solved by calculating the Green's function in NLO. Numerical results for the Large 
Hadron Collider at = 14TeV and 7TeV are presented for various characteristic 
scenarios. 



1 Introduction 



The search for new particles, predicted in supersymmetric models, is one of the important 
tasks of the experiments at the Large Hadron Collider (LHC). The detailed determination 
of particle masses and couplings will be crucial for the discrimination between various 
manifestations of supersymmetry (SUSY) and alternative models, even more so if one 
wants to distinguish between the different variants of supersymmetric models, to identify 
the origin of breaking of supersymmetry and to measure the model parameters. One of 
the SUSY signals will be events with missing energy or missing transverse momentum, 
resulting from cascade decays of squarks and gluinos into the lightest supersymmetric 
particle (LSP) which escapes detection. The existence of squarks, gluinos and the LSP is 
definitely a key prediction of supersymmetry. 

Depending on details of the models, in particular the masses of squarks and the LSP, 
gluinos with masses up to 3 TeV [1] could be detected. The search strategy is quite 
different if gluinos are heavier than squarks or if they are lighter. In the former case 
{rrig > rriq) the two-body decay into a squark and an anti-quark (or its charge conjugate) 
dominates, with the subsequent decay of the squark into a quark plus a chargino or 
neutralino. In the latter case [rag < rriq) the two-body channel is kinematically forbidden 
and the now dominant three-body decay into quark, anti-quark and neutralino or chargino, 
mediated by the virtual squark, leads to a small decay rate. In the limit of extremely 
heavy squarks, a scenario denoted Split SUSY, the gluino is quasistable and hadrons 
composed of gluinos and gluons or quarks may travel macroscopic distances. The search 
strategies will be markedly different in the various cases. 

The importance of squark and gluino searches has motivated a series of detailed studies 
of hadropro duct ion cross sections for squarks and gluinos. The lowest order has been 
evaluated long time ago [2H1]. Subsequently the next-to- leading order (NLO) SUSY-QCD 
corrections were calculated [5H7] , more recently the effect of soft-gluon resummation [8HT3] 
was included. The present paper will be concerned with gluino-pair production close to 
threshold, which exhibits a number of peculiar features. 

As a consequence of their colour-octet representation the production cross section of 
gluinos is large and perturbative corrections are particularly important. Furthermore, 
the threshold region is strongly affected by final state interaction, which in leading order 
is related to Sommerfeld rescattering corrections and which, compared to the similar 
situation in top-anti-top production [15], is amplified by the ratio {Ca/CfY = (9/4)^. 

Corresponding to the different decay modes and rates, two complementary scenarios must 
be considered: For relatively stable gluinos (corresponding to the case rrig < rriq) gluino 
pairs may form non-relativistic boundstates, denoted gluinonia, which decay through 
gluino pair annihilation into a pair of gluon jets or a quark plus an anti-quark jet. This 
possibility has been described originally in Refs. [TBHTU] where the basic features like 
boundstate quantum numbers, leading terms in the potential, spectra and some of the 
production mechanisms were investigated. More recently this aspect has been studied in 
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Refs. [201422] . Relatively stable gluinos are motivated by the proposal of Split SUSY [231 
[21] (see also Refs. [25||26]). which suggests heavy squarks and, correspondingly, long-lived 
gluinos. A detailed study is presented in [27], which gives the higher order corrections to 
the potential, to the boundstate energy spectrum and to the production cross section for 
the boundstate in the colour singlet configuration. 

The second scenario is relevant for gluinos with a larger decay rate, of order one GeV or 
higher. In this case the boundstate decay proceeds through the decay of the constituents, 
and for decay rates of several GeV no well-defined boundstates exist. Even in this case 
final state interaction leads to a significant lowering of the effective production threshold, 
an enhancement of the cross section and a strong distortion of the differential cross section, 
in particular of the distribution in the invariant mass of the gluino, with details depending 
on the gluino mass and decay rate. 

This scenario has many similarities with hadronic top quark production close to threshold 
[2Hll2n]- In particular the distribution in the invariant mass of the gluino pair can be 
treated with similar methods. In Ref. [30] the dominant contributions and the leading 
radiative corrections were classified and evaluated and a sizable enhancement of the cross 
section was observed. However, for a complete NLO evaluation the full ISR and hard 
corrections must be included. In the case of the top system these were available from the 
literature on the production of non-relativistic colour singlet and octet quark-anti-quark 
boundstates [T5l|3T] and could be directly applied to the case of unstable top quarks. 

In contrast, the situation is more involved in the case of gluino pairs. On the one hand the 
production cross section depends in addition on the squark mass; for gluon induced ampli- 
tudes in NLO for quark-induced amplitudes even at tree level. On the other hand a pair 
of gluinos can be combined into boundstates transforming under a variety of irreducible 
representations partly with attractive, partly with repulsive interaction. In Ref. [30] the 
hard correction was extracted from the NLO result for open gluino production. However, 
this continuum cross section corresponds to a weighted sum of the different contributions. 
In the present paper these corrections will be evaluated separately for all the relevant 
representations, together with the individual ISR corrections. This allows to obtain the 
cross section including the full NLO corrections in the threshold region, the central topic 
of this work. 

The paper will be organized as follows: For a self-contained treatment we recall in Section[2] 
the quantum numbers of the boundstates, discuss various SUSY scenarios and present the 
qualitative features of threshold production for the case of interest, i.e. for gluinos with 
decay rates comparable to the level spacing of the would-be boundstates. 

In Section lXn we will present the threshold enhancement (or suppression) for the vari- 
ous colour configurations using Green's functions in NLO approximation. These will be 
evaluated similar to those of the ti system discussed in Refs. [28l|29]. In Section l3T2] the 
NLO pair production cross section, i.e. real and virtual corrections will be derived for 
the different colour and spin configurations. Only S waves will be considered. Issues of 
renormalization, in particular the usage of dimensional reduction (DRED), and the role 
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Table 1: Lowest-lying states of the gluinonium spectrum. L and S correspond to the 
angular momentum and spin quantum numbers, respectively. 



of virtual gluino or squark pairs will be discussed, together with the choice of the proper 
value of the strong coupling a^. 

Using this input, the hadronic production cross section can be evaluated in a straightfor- 
ward way in SectionHl We will limit the discussion to proton-proton collisions at 7 and 
14 TeV and give results for several of the SUSY scenarios discussed in Section[2l Section|5] 
contains our conclusions. Appendix |X] contains useful relations for the generalized hyper- 
geometric function 4F3, and some of the longer formulae are relegated to Appendix [Bl In 
Appendix O details of the benchmark scenarios used in this paper are provided. 

2 SUSY scenarios, gluino boundstates and threshold 
behaviour 

Let us briefly recall the quantum numbers of gluino pairs in the threshold region, classified 
according to their colour, spin and orbital momentum configurations [T614T8] . They differ 
from those of quark-anti-quark states due to the restrictions arising from the Majorana 
nature of gluinos, and due to their different colour assignment. 

Two colour-octet states can be combined into irreducible representations as follows (see 

e.g. m) 

8® 8 = 1, ©8, ©8a©10a©T0a©27,, (1) 

where the subscript indicates (anti-)symmetry with respect to their colour index. Fermi 
statistics and the Majorana nature of the gluinos lead to additional restrictions. For the 
symmetric colour configurations 1^, 8^ and 27^ antisymmetric spin-angular momentum 
wave functions, (—1)'^+'^ = 1, are required, for the antisymmetric colour configurations 
8a, IOq and IOq symmetric ones, (—1)-^"'''^ = —1. The intrinsic parity of a Majorana particle 
can be chosen to be imaginary leading to negative intrinsic parity of the boundstate. For 
a few lowest orbital angular momenta the boundstate quantum numbers are listed in 
Tab. [TJ The transformation of these states under charge conjugation is more involved^ 

^We thank Y. Kats and D. Kahawala for drawing our attcntfon to this issue. 
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Since we restrict the discussion to boundstates, only the colour configurations 1^, 8s and 
8a with attractive potentials are discussed in the following. Colour quantum numbers now 
play a non-trivial role. A self-conjugate Majorana particle (without colour) transforms 
with C = +1 under charge conjugation, the transformation of coloured constituents and 
their boundstates, however, depends on their colour index. Let us consider boundstates 
of two Majorana particles. For the symmetric singlet states with the colour projector 6ab 
the charge parity of both constituents is identical, leading to overall positive charge con- 
jugation. This is consistent with the fact that the decay of the pseudoscalar boundstate 
into two photons is non- vanishing [27j, and agrees with earlier statements in the litera- 
ture [T6] - [T8] . The quantum numbers of the antisymmetric octet ground state (^Si) with 
= 1~ can be obtained by considering the wave function of the boundstate ~ Qbl'^Qcfahc 
(where cjb and cjb are the Majorana fields) with C-quantum numbers identical to those of 
the gluon. These, in turn, can be obtained from the relation [33l|3l] between the gluon 
field and its charge conjugate A'^^ 

V^Al = -TlAf, (2) 

where Ta are generators of SU(3) in a specific representation, the corresponding trans- 
posed matrices. Specifically this implies C = —1 for with a = 1, 3, 4, 6, 8 and C = +1 
for a = 2, 5, 7 and correspondingljH for the boundstates with = 1~. The same assign- 
ment is also valid for the {gg)a state with = 

The charge parity of the pseudoscalar state {gg)s in the symmetric octet representation 
which corresponds to the wave function ~ gbl^gcdabc is given by C = +1 for a = 1, 3, 4, 6, 8 
and C = —1 for a = 2, 5, 7 and hence opposite to the one of the antisymmetric statesjf] 
Note that for this assignment C-parity is conserved in the decay of the boundstate to two 
gluons which proceeds through the coupling ~ dabcA^A'^. The same assignment is valid also 
for the other states in the symmetric octet representation. In total this can be summarized 
by defining a factor Ca with Ca = +1 for a = 1,3,4,6,8, Ca = —1 for a = 2,5,7, and 
assigning charge conjugation —Ca for the gluon field charge conjugation +Ca to the 
gluino, —Ca to the antisymmetric octet and +Ca to the symmetric octet boundstate. 
Restricting ourselves now to the near threshold region, only 5'-wave configurations will be 
retained. 

Depending on the representation R of the boundstate, the interaction between the two 
gluinos can be either attractive, repulsive or absent (in lowest order). In lowest order the 
coefficient C'^^ of the QCD potential which governs the final state interaction is given by 
the expectation value of the product of the colour generators F^-F^i, taken between two- 
particle states in the respective representation. This product, in turn, can be expressed 
by the eigenvalues of the quadratic Casimir operator of the constituents, Ca = 3, and the 

^This result is at variance with [l8l[27] which find C = +1 for all colour labels a as well as with [35ll36] 
which find C = — 1 for all values of a. 

■^This result is at variance with [TS } [27 1 135 1 136 ] which find C = +1 for all colour labels a. 
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Table 2: Colour interaction of two SU(3) octets, 
boundstate in representation R, Cr = [F^^: 

^[R] ^ pa,l . pa,2 ^ }_ ^^pR^ _ ^pa,iy _ (^a,2)2j =^{Cr- 1C a) ■ (3) 

The results are listed in Tab. [2J For the cases with negative (positive) coefficients, corre- 
sponding to attraction (repulsion), the cross section will be enhanced (suppressed). Also 
the NLO correction of the QCD potential, which will be needed below, is proportional to 
the same group theoretical factor [30] . The classification described in Tabs. [1] and [2] is of 
course also applicable to continuum production and will be important for the description 
of final state interaction. 

As mentioned in the Introduction and discussed in the literature [T6HT81 127 1 130] the phe- 
nomenology of gluino pair production in the threshold region is governed by the relative 
size of the decay rate of a gluino, Tg, compared to the rate for gluinonium annihilation 
into gluons, Tgg, on the one hand, and by the relative size of Tg and the level spacing AM 
between the ground state and the ffist radial excitation of the colour singlet boundstate on 
the other hand. Specifically we adopt for this comparison the single-gluino decay rate of 
the boundstate, corresponding to 2Tg. The choice of AM = \Ei — is motivated by the 
fact that the binding energy per se depends evidently on the choice of the mass definition 
(pole mass, potential subtracted mass, . . . ) while AM is convention independent. 

If the decay rate 2rg is smaller than the annihilation rate of the (S'-wave) bound- 
state, which in lowest order is approximately given by Tgg ^ C^a^ |_R(0)|^ /(2m?) ^ 
(Cyias)^mg/4, the signatures of boundstate and open gluino production are distinctively 
different: Boundstates produced below the pair threshold decay into two gluon jets (no 
missing energy), in contrast to the complicated cascade decays of gluinos above threshold. 
This case (class A), evidently true in Split SUSY, was discussed in detail in Ref. [27], with 
emphasis on the production of colour singlet boundstates. 

If Tgg <, 2rg, tliB coustitueuts decay before they can annihilate and the qualitative decay 
signatures of bound and unbound gluinos are practically indistinguishable. Nevertheless, 
the distribution in the invariant mass of the gluino pair will still be modulated and strongly 
affected by final state interaction. As long as 2rg < AM, this mass distribution will still 
reflect the presence of the boundstates, at least the enhancement resulting from the lowest 
IS* resonance with quantum numbers O""*" (class B). For even larger gluino decay rates 
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i.e. 2r^ > AM these structures have essentially disappeared resulting in an unstructured 
threshold behaviour (class C). Nevertheless the cross section is still modified by final 
state interaction and receives a significant contribution from the region below threshold 
as far as colour singlet and octet rates are concerned and a sizable suppression for the 27 
representation. 

The dependence of Tgg, 2Tg and AM on nig is illustrated in Fig. [H In Fig. [1^ we compare 
AM and Tgg, evaluated in NLO [27], to the decay rate 2Tg, evaluated for three generic 
choices of squark masses 0.5 TeV, 1.0 TeV and 1.5 TeV. In Fig. [T|d fixed ratios between 
rrig and rriq {rrig/mq = 0.75, 0.90, 1.05, 1.25) are adopted. In Tab. |5] (Appendix [C|) we 
define 17 benchmark points on the basis of various SPS scenarios introduced in Ref. [37] . 
Decay rates and SUSY masses corresponding to our benchmark points are calculated 
with the programs SuSpect [SB] and SDECAY [22] • The essential information, i.e. gluino 
mass, averaged squark mass, level spacing of gluinonium, gluino decay rate, gluinonium 
annihilation decay and dominant gluino decay channels, is listed in Tab. [31 Tabs. |6]and[7] 
in Appendix [C] contain detailed information on the squark masses of the benchmark 
points. Gluino masses and decay rates, level spacing and annihilation rate of gluinonium 
for these benchmark points are shown in Fig. [lb. It is clear from this figure and Tab. [3l 
that all three possibilities, corresponding to class A, B and C could arise and should be 
discussed. Note, that even the three cases which fall into class C are fairly close to the 
boundary between B and C, such that a structured threshold behaviour is typical for all 
benchmark points discussed in this paper. Class A has been studied in detail in [27], with 
emphasis on colour singlet production. The present paper will be concerned with class B, 
with results also applicable to class C. 

From these considerations, from Fig.l and from the discussion presented below it is clear 
that the parameter p = mg/rriq which characterizes the relative size of gluino vs. squark 
mass and thus the gluino decay rate is decisive for the assignment to class A, B, or C. 
A value close to one is characteristic for models at the borderline between A and B, a 
value around 1.14 corresponds to models at the borderline between B and C. Benchmark 
points a and q with p = 1.3, 1.12 and 1.16 respectively, although excluded by recent 
LHC results [M12] alreadjfl nevertheless provide important insight into the structure 
of the threshold behaviour and the relative importance of the different spin and colour 
configurations. This will be illustrated by two scenarios X and Y with the parameters 
listed in Tabs. 3 and 5 and which are not excluded by LHC results. Model Y corresponds 
to point 10.3.3 from Ref. [12], in model X the parameters mo and mi/2 have been chosen 
such as to ensure a relatively small gluino decay rate. The phenomenological discussion of 
gluino pair production for these two scenarios will be presented at the end of Section 4. We 
restrict ourselves to configurations suggested by the Constraint Minimal Supersymmetric 
Standard Model, the results, however, depend mainly on the gluino decay rate only and 
thus are valid in a wider context. 

"'For the Constraint Minimal Supersymmetric Standard Model limits of around 1 TeV are quoted for 
rrig and m^. Note, however, that the calculation performed in this paper is more general and can be 
applied to less restrictive supersymmetric scenarios. 
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Figure 1: Level spacing AM and annihilation rates Tgg (solid curves) compared to the 
gluino decay rates 2rg. Fig. (a): for three different squark masses (dashed: 0.5 TeV, 
dash-dotted: 1 TeV, dotted: 1.5 TeV). Fig. (b): for fixed ratios (p = mg/rriq) between 
squark and gluino masses. Fig. (c): for the benchmark points (see Tab. [3]). 
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benchmark 
point 


mg [GeV] 


mg [GeV] 


AM [GeV] 


2rg [GeV] 


r,, [GeV] 


class 


dominant decay 
channels 


(a) 


606.11 


541.04 


13.89 


9.08 


0.37 


B 


g bj), bib 


(b) 
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12.04 


6.28 
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B 


g C2C, C2C, 

"5,2 U, 
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g — > bib, bib 
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0.40 
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g —5- bib, bib 


(e) 


826.71 
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0.43 
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Table 3: Comparison of gluino masses, average squark masses, gluinonium level spacing, 
single decay rate, annihilation rate and dominant decay channels for the 19 benchmark 
points defined in Appendix [Cl Tab. [5l 
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Figure 2: Feynman diagrams contributing at LO to gg — ?■ gg and qq — gg. 



Let us now recall the qualitative aspects of the production mechanism. Similar to the 
case of top quark production [28l[29] (see also [13]) the cross section for a bound state T, 
differential in M, the invariant mass of the gluino pair, can be decomposed into a factor 
representing the hard, short distance part of the production process, and a factor given 
by the imaginary part of the Green's function evaluated at the origin and the convolution 
with the luminosity functions 
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As usual s and S" denote the partonic and the hadronic center-of-mass energy squared, 
respectively, and r = s/S. The lower limit of the r integration is given by p = M'^/S. 
The superscript of the Green's function refers to the colour state of T and fip and fi^ 
denote the factorization and renormalization scale. 

In leading order (LO) gluino pairs can be produced by gluon fusion or quark- ant i- quark 
annihilation (Fig. [2]), at NLO also the quark-gluon channel contributes. 

The strength of final state interaction at threshold is determined by the respective colour 
representation of the gluino pair. To disentangle the different representations it is conve- 
nient to use the following projectors 
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where Pio is the projector on 10a © 10a and 



SacSM^x'"' = Nx, 
Kcd 



Sac^bd ■ (7) 



X 



The normahzation constant Nx is given by the dimension of the representation X, where 
TVio = A^'io^ + NjQ^ = 20. As stated above, colour representation, spin and orbital angular 
momentum of the gluino pair are interlocked. S-wave gluino pairs in symmetric and anti- 
symmetric colour configurations combine into spin singlet and triplet states, respectively. 
For the quantities J^ij^r (as introduced in Eq. ([5])) describing the hard kernel one finds 
in Born approximation 

^1^. = (8) 
with the non-vanishing normalization factors 

ATlJ] = 1,2,3 for T=l„8„27,, 

Here z = M'^/s and r = m^/m'- where rfiq is the averaged squark mass as introduced in 
Tab. [3l The Green's function which depends on the energy E = M — 2mg and the decay 
rate of the gluino Tg is obtained from the Schrodinger equation 



«V)^ , [1,8,10,27]/^. 

+ Vc 



m-g 



{E + tTg)\G^''''''^^^'\f,E + tTg) = (5(=^)(f) (10) 



with the "Coulomb" potential V^r (f) = -C^^^as{i2R)/r (and Cl^l 
3,3/2,0,-1 for R = 1,8,10,27) in lowest order. NLO corrections to the hard 
kernel, the potential and the Green's function will be discussed in the next section, 
numerical results will be presented for a variety of scenarios in Section |H From now 
on we shall limit our discussion to three typical benchmark points (p), (a) and (q) (see 
Fig. [Hand Tab. [3]). For (p), due to the close proximity of squark and gluino masses, the 
decay rate of the boundstate, which is approximately ^Tg, is significantly lower than AM 
and, as a consequence, the enhancement from the lowest lying resonances is well visible. 
This is the case we will discuss in most detail. The other two points serve to illustrate 
the case of gluinos with somewhat larger decay rates. 

The imaginary parts of the LO and NLO Green's functions (for R = 1,8, 10, 27) are 
displayed in Fig. |3]^a). We have adopted twice the Bohr radius as characteristic scale 
in as for the attractive potential, as{Mz) = 0.1202 as starting value and two-loop 
running. This leads to ai^^ = a,(C[^]aPWg) = {0.1034,0.1124} for [R] = 1,8 and 
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Figure 3: Imaginary part of the Green's functions for benchmark points (p), (a) and (q) for 
the singlet (dashed), octet (dash-dotted), decuplet (sohd) and twenty-seven configuration 
(dotted). For the decuplet we show the free Green's function. LO and NLO curves for 
twenty-seven lie on top of each other. 
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Figure 4: LO prediction for the invariant mass distribution for scenario (p). 



as = cisi 

\cm 

of rrig) = 0.1184 for the repulsive potential of the twenty-seven rep- 
resentation. Two further characteristic examples are shown in Figs. [3]^b) and (c) cor- 
responding to benchmark points (a) and (q) with significantly larger F^. LO and NLO 
predictions are shown for the singlet (dashed) and for the octet configuration (dash- 
dotted), anticipating the results of Section 13. 1[ In our approximation the potential of 
the decuplet vanishes and thus the free Green's function is shown (solid). For the repul- 
sive twenty-seven configuration the NLO corrections are quite small, i.e. LO and NLO 
nearly coincide. The accidental near degeneracy of the first radial excitation from G^^' 
with the lowest enhancement of G^^' follows trivially from the ratio between the strengths 
of potential (C'^^/Ct^')^ = 1/4 and the excitation spectrum for the Coulomb potential. 
En ~ l/n^. For (p) the lowest resonances are still nicely separated, for (q), which is 
close to the border between class B and C (see Fig. [11(c)), individual resonances have 
nearly disappeared. Nevertheless final state interaction leads to a significant modifica- 
tion of the threshold behaviour and to marked differences between the different colour 
representations. The Green's functions for symmetric and antisymmetric octet are obvi- 
ously identical and denoted by G^^\ the difference between the two states is only the spin 
configuration. 
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The LO prediction for the cross section is shown in Fig. H] separately for the four different 
contributions and the sum. We have employed the PDF set MSTW2008LO [15] which 
fixes the value as{Mz) = 0.1394. We use two-loop running and one-loop decoupling (as 
implemented in the program RunDec [16]) to obtain the strong MS coupling at the scale 
fi = 2mg. In a next step we transform the coupling to the DR scheme in the full SUSY 
QCD theory with one-loop approximation (see, e.g., Ref. [IZ]). Note, that the qq induced 
production of 8^ is strongly suppressed as a consequence of the destructive interference 
between the amplitudes with s-channel gluon and t- as well as u-channel squark exchange 
and the near degeneracy of squark and gluino masses 

2 2\ 2 

-i 1 ^ 3.5- 10"^6.1■ 10"^l.l■ lO^S (11) 

m| + m| / 

for the benchmark points (p), (a), (q). Close to threshold the repulsive final state inter- 
action of 27s leads to a strong suppression of the cross section. For larger energies this 
effect disappears quickly as a consequence of the large multiplicity of states in the 27^ 
representation. 
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3 Next-to-leading order corrections 

The NLO corrections to the cross section, whose evaluation is the main subject of this 
work, can be separated into those for the Green's function Im {G(0, E + iVg)} and those 
for the hard coefficients J^ij^T-, as defined in Eq. ([5]). 



3.1 Green's function 



Following the idea of the Green's function formalism developed in Refs. [181 - I50] we start 
with using the NLO potential in momentum spac^ 



V!?' {q \ 
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^[27] 


= -ic 
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3 
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(13) 



"Tor colour triplet or octet states combined into singlet boundstatcs the potential is even known to 
NNNLO [511l54| . for the combination of colour triplet states into octet representations see |55| . 
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and 



31^ 20 



(14) 



We have checked by an exphcit calculation that for all colour configurations the same 
coefficient ai is obtained. The renormalization scale relevant for potential and Green's 
function is denoted by fic and chosen to be the solution of fic = \C^^^ \ as{fJ,G)^g- The 
strong coupling is defined in the MS scheme and for five active fiavours. In principle 
one might include top-quark mass effects, however, at the present level of precision these 
terms are still irrelevant and are neglected in the potential. 

The Green's function for the top-anti-top system is known in compact analytic form [56] 
(see also [57]) and the result is easily applied to the present case 
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The n-th derivative ip^"^^ = ip^^^l — n) of the digamma function ijj{z) = 'jE + id/dz) lnr(2;) 
is evaluated at {1 — k) and the Generalized Hypergeometric Function 4F3 is defined in 
Appendix Rl 

Solving the Schrodinger equation perturbatively induces poles in the Green's function of 
the form ja^E'^^/ [E^^ — {E + iTg)'j } which become large in the vicinity of E = E^'^ 
for a small decay width Tg. To obtain a proper Green's function with single poles one has 
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to resum the muhiple poles as proposed in [56] by adding the term 

F^^ (1 + g JO 
£;LO(i + eia,)-(E + zr^) 

r ^Lo 

which is of order a^. The Schrodinger wave function at the origin as well as the binding 
energy are given by \^n{0)\'^ = F^^ (1 + aji + . . . ) and E„ = (1 + a.d +...)• 
The quantities F^^, /i and ei can be found in Appendix A of Ref . [27] . 

Attention has to be paid to the numerical evaluation of the Generalized Hypergeometric 
Function 4F3 in Eq. f|T5|) . This function has a branch cut for its last argument on the 
real positive axis starting from 1. Hence the series defining this function converges for 
Re (1 — k) > 1, a condition potentially violateci§ for small \E + iTg\. In Appendix|X] it is 
shown how to circumvent this problem by a suitable transformation^ 

Due to the factorization, formulated in Eq. (jl]), the renormalization scales of Green's 
function fie and short distance part /i^, appearing below, can be chosen independently. 
For our choice of /ig, twice the inverse of the Bohr radius of the two gluino system, 
the Coulomb- Green's function has a well-convergent perturbative series as discussed in 
Ref. [58] for the top-anti-top pair in the colour-singlet configuration. 



+ 



{Elf^-{E + iVg)f Elf^-{E + iTg) 



3.2 Short distance corrections 

The calculation of the partonic cross sections has been performed in dimensional regular- 
ization (DREG). The Feynman diagrams have been generated using FeynArts [59l|60], for 
the evaluation we have used in-house FORM [61] programs and FormCalc [60]. Masses 
and wave functions are renormalized on-shell while minimal subtraction (MS) is used 
for the coupling. DREG violates supersymmetry because of a mismatch of the degrees of 
freedom of the gauge bosons and their fermionic superpartners. This problem has been re- 
solved by introducing dimensional reduction (DRED) [62] as regularization method where 
the elegant features of DREG are maintained. While the on-shell masses do not depend 
on the choice of scheme the coupling has now to be calculated in DRED combined with 
minimal subtraction, the so-called DR scheme. The transition between the MS and DR 
parameters is known up to two loops [631IM] . For our calculation the one-loop relations 
are sufficient which we need for the Yukawa coupling, which appears in the gluino-quark- 
squark vertices, and the gauge coupling in the three-gluon, gluon-quark and gluon-gluino 
vertices. Our final results are expressed in terms of = ai^^*"^'', the strong SUSY QCD 
coupling in the DR scheme. 

6 The condition is fulfilled for 1 > \e];'^ \ {\/WTf^ - E) / [2{E'^ +T'^)\. Hence this problem does 
not arise for the tt system. 

^Wc thank Yuichiro Kiyo for communications concerning this point. 



16 



The corresponding amplitudes have to be projected onto the proper spin configuration 
[65|l66] and the corresponding colour representation (see Eq. (|6])). As discussed in Section [2] 
the colour-symmetric configurations (with L = 0) have spin 5 = 0, the antisymmetric 
ones S = 1. For the processes, which are non- vanishing in Born approximation, the hard 
part of the partonic cross sections can be written in the form (with the factors J\flp given 
in Eq. dU) 
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The remaining ones are given by 

ij-t-T 
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Furthermore, J^gq^r is equal to Tgq-^T- The quantities Vj^ and T^j^- denote the virtual 
and real corrections and are obtained from the results listed below (see Eqs. (pT]) - ( 125|) ) 
by dropping the infrared singularities, which cancel between \>\- and 7?. . For the Born 
approximation in ci = 4 — 25 dimensions a factor (1 — e){\ — 2e) for the pseudoscalar and 
(1 — e) for the vector states has been taken into account. (Note that collinear singularities 
from T^^j are already absorbed in the PDFs.) 

The only non-vanishing contributions to the hard virtual corrections read 
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where mt = 172 GeV, /.(Q^) = {An Q^Y T {1 + e) and /3o^^° = 3Ca - ^TpUf (with 
72/ = 6) is the one-loop coefficient of the SUSY QCD beta function and the scalar functions 
hi, b[ and C5 are defined in Appendix [Bl The limit mt ^ has been taken wherever 
possible. 

The gg initiated processes into 1^, 8^, 27^ receive contributions from the diagrams depicted 
in Fig. [5l which do not involve squarks and from diagrams with virtual squarks depicted 
in Fig. El The latter are proportional to Uf, the number of quark and squark flavours. For 
the squarks we assume equal masses. (The general calculation for different masses for left- 
and right-handed squarks and for different generations is straightforward, however, the 
formulae are lengthy. Explicit results are presented in Ref. |67].) The qq initiated process 
into 8a receives contributions from diagrams without and with squarks in LO already (see 
Fig. [2]). Consequently also the corrections contain contributions from diagrams without 
squarks (independent of r) (Fig. [7]) and with squarks (Fig. [S]). All squark mass dependent 
terms are collected in the function A^^^^ which is defined in Appendix iBl 

Born term and virtual corrections are absent for decuplet production in gluon fusion 
and quark- anti-quark annihilation (see Eq. ([H])). The corrections from real radiation as 
discussed below give a (small) non- vanishing result. 

IT] 

The quantities TZ\j result from the real corrections to the processes and contain the 
subtraction terms which originate from the renormalization of the PDFs in the MS scheme. 
Thus the remaining infra-red divergences exactly cancel the ones present in the virtual 
corrections of Eq. f l2T|) . Let us start with the process gg — )■ Tg which does not involve 
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Figure 5: NLO contributions to the gg initiated processes without virtual squarks (dotted 
hues correspond to ghost fields). 
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Figure 8: NLO contributions to the qq initiated processes with virtual squarks. 
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squarks (see Fig. [9]). For the individual colour configurations we obtain 
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where nj = 6. The functions Vij are listed at the end of this Section. The conventional 
plus-distributioi]0 is employed to regularize the singularity at 2 = 1. Note, that the cross 



^The plus-distribution follows the prescription dz 
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Figure 10: NLO contributions from gq initiated processes. 



section for the decuplet configuration (plus a gluon) is non-vanishing, albeit small. 

The results for the gq and gq initiated processes gq — )■ Tq and gq — Tq are identical and 

IT] IT] 

we only list the former. The corresponding function TZgg can be split into a term Tigg , 
which originates from the diagrams depicted in Fig. [Tor a)-(c). and is independent of the 
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Figure 11: NLO contributions to the qq initiated processes. 



squark mass and a second term J^]^' from the diagrams depicted in Fig. [TOl(d)-(h) and 
their interference with those of (a)-(c). The dependence on the factorization scale is 



contained in . One finds 
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r-i 1 [■'T-ti 

and T-Lgq of course coincides with the result of Eq. (12) of Ref. |2^. The functions J-'gq 
are listed in Appendix El see Eqs. ( I33|) - (l371) . A similar decomposition can be made for 



25 



the qq initiated processes (cf. Fig. [TT 
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and J-"'^' is again listed in Appendix [Bl The sphtting functions entering the previous 



expressions are defined as 
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4 Hadronic production 

As described in Eq. (j4]), the luminosity function Eq. ([5]) is convoluted with the partonic 
cross section, which in turn is composed of the Green's function (Eq. ( !T5|) ) and the short 
distance corrections (Eqs. ( IT9l) and (|20l) ). For the numerical evaluation we use the PDF 
set MSTW2008NLO [IH] which corresponds to as{Mz) = 0.1202. (Note, that for the 
LO study a value as{Mz) = 0.1394 was adopted.) From this value as starting point we 
compute a^^^^ employing the same procedure as described before Eq. ffTTj) . We thus 
arrive at the cross section for proton-proton collisions, which is still differential in the 
mass of the gluino pair. In a first step we restrict ourselves to scenario (p) and display 
the region very close to threshold (Fig. [T^ a)) and a wider range with M^^^ = 2000 GeV 
corresponding to f^^^^^ = 0.72 (Fig. [T2r b)). The contributions from the different colour 
configurations are shown separately. Close to threshold the strong enhancement of the 
singlet is contrasted with the strong suppression of the twenty-seven representation. For 
larger M and correspondingly larger relative velocities the final state interaction becomes 
less relevant and about 200 GeV above threshold the twenty-seven representation starts 
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Figure 12: NLO prediction for the differential cross section for scenario (p) for two different 
regions of M. 
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Figure 13: Prediction for the differential cross section in NLO using the Green's function, 
in comparison with the fixed order cross section without and with vanishing single decay 
width for scenario (p). 



to dominate. (This behaviour is quite similar to the one of top pair production with 
dominant singlet close to and dominant octet far above threshold.). The process qq — )■ [8a] 
is non-vanishing in Born approximation. However, for squark and gluino masses being 
roughly equal it is strongly suppressed. This suppression is also present at NLO. The 
contribution from the decuplet, which is absent in Born approximation, remains small 
throughout. In Fig. [13] we compare this result (solid curve) with the NLO result using a 
fixed order treatment for the Green's function, i.e. we replace the imaginary part of the 
Green's function by its expansion in a^, first keeping F^ non- vanishing (dotted curve) and 
then in the limit Fp — )■ (dashed curve) 

l^G^R] ^ fi + ct^l^) , (27) 

leaving hard correction and PDFs unchanged. In the threshold region the latter choice is a 
valid approximation to the complete fixed order NLO prediction for the cross section. For 
invariant masses above 1700 GeV fixed order and Green's function modulated approach 
agree reasonably well, between threshold and 1600 GeV they differ significantly. (Inclusion 
of finite width effects is quite irrelevant in the fixed order treatment.) The integrated 
difference between dashed and solid curves amounts to 0.183 pb and is a measure of the 
threshold enhancement, that would escape the strict fixed order treatment. Relative to 
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NLO calculation 
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approximation 


[fb/GeV'] 


0.0954 


4.29 
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12.9 


[fb/GeV^] 


0.0862 


6.50 


14.0 


0.136 
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gNLO/^LO _ 1 




51.5 


63.4 


16.2 




82.2 



Table 4: Comparison of full and approximate NLO results (see text). 



the total cross sectioij^ for hadro-pro duct ion of two gluinos, which, within scenario (p), 
amounts to cTtot = 2.59 pb, this corresponds to an enhancement of 7.1%. 

The renormalization {^r) and factorization {^f) scale dependence of the differential cross 
section is shown in Fig. [Tl]for three different choices, namely = = ^gi'^^gi'^f^g 
and again for the two regions in M and for the LO and the NLO prediction. Note, that 
we keep the Green's function unchanged and identify factorization and renormalization 
scales. For the LO prediction we use the parameters described at the end of Section 2. 
Between LO and NLO we observe a slight shift of the location of the resonance peaks by 
about 10 GeV, a slight increase of the cross section by about 15% and a reduction of the 
scale dependence. For the NLO prediction this residual /i dependence amounts to 15%. 

In Figs. [15] and [16] the corresponding results are shown for scenarios (a) and (q). The 
gluino masses of (p), (a) and (q) are comparable, the difference between gluino and squark 
mass, however, increases and, correspondingly, the gluino decay rate. For scenario (a) with 
AM = 13.89 GeV and 2rg = 9.08 GeV the 15 peak is still clearly visible, for scenario (q), 
with AM = 15.79 GeV and 2r^ = 22.92 GeV, the resonant structures have essentially 
disappeared. Nevertheless, final state interaction leads to a significantly enhanced cross 
section in the threshold region also in these two cases (Fig. [TTT a) and (b)). 

The results for the invariant mass distribution with ^/s = 7 TeV are shown in Fig. [TST a) 
and (b), restricting ourselves again to scenario (p). These results are qualitatively similar 
to those of Fig. [12] A reduction of the cross section by a factor 20 is observed. 

In Tab. [4] we investigate the relative size of the corrections for the different states and 
contrast our result with those based on an approximation discussed in [^iOj. The latter 
employs a hard correction factor, which is different for quark-anti-quark annihilation 
and gluon fusion, does not distinguish between different colour representations and has 
been numerically extracted from the continuum result [5]. Furthermore, it includes the 

^This result has been obtained with the code Prospino2 [5j in NLO, using CTEQ5 PDFs. 
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Figure 14: Renormalization- and factorization-scale dependence of the differential cross 
section for scenario (p) for two different regions of M. 
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Figure 15: Renormalization- and factorization-scale dependence of the differential cross 
section of scenario (a) for two different regions of M. 
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Figure 16: Renormalization- and factorization-scale dependence of the differential cross 
sections of scenario (q) for two different regions of M. 



32 



0.015 



0.012 



> 

^ 0.009 



D 



0.006 



0.003 





1 1 i 1 1 i 1 1 1 1 1 i 1 


1 1 1 1 1 1 1 

Green . 






FO ( Tg = 4.54 GeV ) 


- 




FO(rg = 0GeV) - - 




/ ^ 

/ 

/ 






1 

i 

i 




/ 


■i (a) 

mj = 606.11 GeV 

• 1 8 

/ rriq = 541 .04 GeV 

' 1 2rg = 9.08 GeV 
1 ... 1 ... 1 ... 1 





1200 



1300 



(b) 



1400 1500 

M [GeV] 



1600 



1700 



> 



c 

"D 



0.004 



0.003 



0.002 



0.001 



Green 

FO ( Tg = 1 1 .46 GeV ) 
FO ( Tg = GeV ) 




(q) 

mg = 719.66 GeV 
mp = 618.86 GeV 
2rs = 22.92 GeV 



1400 



1500 



1600 



1700 

M [GeV] 



1800 



1900 



2000 



Figure 17: Prediction for the differential cross section in NLO using tlie Green's function 
in comparison witli the fixed order cross section without and with vanishing single decay 
width for scenarios (a) and (q). 
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Figure 18: NLO prediction for the differential cross section at y/s = 7TeV for scenario 
(p) and comparison to the fixed order cross section with vanishing single decay width. 
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logarithmically enhanced terms from soft radiation, which are determined by the colour 
charge of initial and final states. To allow for a consistent comparison, we employ for both 
cases MSTW2008LO PDFs with a,(Mz) = 0.1394 in LO and MSTW2008NLO PDFs with 
as{Mz) = 0.1202 for the (approximate) NLO result. However, following |30] we use for 
the approximate treatment two loop running with n/ = 5 both for LO and NLO. The 
corresponding values are listed in the table. In Tab. IHwe also show the results for the 
quantities and s^^° which were obtained from the corresponding results for da/dM 
with G = mi/{4:T[). In this way we remove the dependence on the Green's function. The 
hard corrections are only weakly dependent on M and we adopt M = 2mg as reference 
point. Furthermore, we adopt benchmark point (a), corresponding to rrig = 606.11 GeV 
and rrig = 541.04 GeV. 

The difference between the two LO results can be traced to the different values of a^. The 
difference between the relative size of the corrections is due to additional subprocesses 
present in the full calculation and the appearance of additional contributions with virtual 
squarks. 

As stated in Section 2, the three scenarios (p), (a) and (q) serve to illustrate the change 
of the excitation curve by moving from relatively smaller gluino decay rate up to a situ- 
ation with 2r^ even somewhat larger than the excitation energy AM. Although already 
excluded by recent LHC results, their detailed discussion serves to illustrate the impact of 
NLO corrections and threshold enhancement. From the phenomenological side, however, 
it seems appropriate to also present results for scenarios still consistent with the LHC 
limits. Indeed, the behaviour is quite similar for the two scenarios (X) and (Y) (see Tabs. 
3 and 5). In view of the smallness of the respective cross section for ^/s = 7TeV only the 
results for a/s = 14TeV are presented. We separate the cross section according the the 
different colour and spin configurations (Figs. [19] and |20|) and compare the results with 
the corresponding ones based on NLO fixed order calculations with vanishing and non- 
vanishing gluino decay rates (Figs. [21] and |22]) . As anticipated, the qualitative features 
of the results are quite similar to those for scenarios (p), (a) and (q). 

5 Conclusions 

The next-to-leading order analysis for hadronic production of gluino pairs close to thres- 
hold has been presented. The matching coefficients were evaluated separately for the 
different colour configurations of gluino pairs with relative angular momentum zero. The 
cross section is strongly affected by final state interaction which is encoded in the NLO 
Green's function and which depends on the gluino decay rate, the colour configuration 
and the invariant mass of the pair. In a first step we have investigated three different 
SUSY-scenarios covering a wide range of gluino and squark masses, and studied the 
renormalization and factorization scale dependence of the result. Compared to the leading 
order prediction the result is more stable, and, for fi between rrig and Arrig, varies by 
±15%. The complete NLO threshold production of (gg) boundstates considered in this 
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Figure 19: NLO prediction for the differential cross section for scenario (X) at 
14 TeV for two different regions of M. 
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Figure 20: NLO prediction for the differential cross section for scenario (Y) at 
14 TeV for two different regions of M. 
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Figure 21: Prediction for the differential cross section in NLO using the Green's function, 
in comparison with the fixed order cross section without and with vanishing single decay 
width for scenario (X) at ^/s = 14 TeV. 



paper enhances the fixed order prediction by typically 7% to 9%. To accomondate the 
most recent limits on squark and gluino masses, we also presented predictions for two 
additional scenarios with squark and gluino masses above 1 TeV. The qualitative features 
of our results remain unchanged. 

A The Generalized Hypergeometric Function 

The Generalized Hypergeometric Function (GHF) is defined by the series 

{bl)n ■ ■ ■ {bq)n n\ 

with the Pochhammer symbols = r(x + n)/r{x) and the restriction 6j 7^ 0, —1, . . . 
for i = 1,2, ... ,q. The series converges if one of the following conditions holds 

(1) p<q,\z\ < 00, 

(2) p = q + l,\z\<l, 

{3) p = q + l,\z\ = l. Re (ELi bn - Et\ > , 
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Figure 22: Prediction for the differential cross section in NLO using the Green's function, 
in comparison with the fixed order cross section without and with vanishing single decay 
width for scenario (Y) at ^/s = 14 TeV. 



(A) p = q+l,\z\ = l,z^l,-l<Re (ELi - E'li «n) < . 



For the numerical evaluation of the Green's function of Eq. f[T5|) it is necessary to evaluate 
the GHF 4-^3(1, 1, 1, 1; 2, 2, x; 1) for a; G C, hence the remaining arguments have to fulfill 
condition (3). For Re(a;) < the convergence of the series is not guarantied. 

In Ref. [68] the following algorithm is introduced which allows to decompose this function 
via partial fractioning in GHFs with convergent series representation. Following Ref. 
we employ the following identity 



4^3(1, 1,1,1; 2, 2, x;!) 



1 



4x^(2 - x) 



4(x -1)^^3(1,1, 1,1; 2, 2, x + l;l) 

+2a;(7-4a;) 4^3(1, 1, 1, 1; 3, 2, x; 1) 
+x(x-2)4F3(l,l,l,l;3,3,x;l) 



(29) 
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For arbitrary values Re(x) < with x ^ —1, —2, . . . this relation is apphed repeatedly, 
leading to the following results 

4^3(1, 1, 1, 1; a, a,x; 1) = 



X 



o?x{x-2{2-a)){a-xY 

a^i^x — lY 4-^3(1) 1, 1, 1; a, a, X + 1; 1) 
+a(a - lfx{'ia + 1 - 4x) 4i^3(l, 1, 1, 1; a + 1, a, x] 1) 
+ (a - lYx{x - a) 4-^3(1, 1, 1, 1; a + 1, a + 1, x] 1) , 



4^3(1, 1,1,1; a, &,a;;l) 



1 



a + 6 + X — 4 



[a-iy 



a{a — b){a — x) 

h{h — a){h — x) 
I - If 

x{x — a){x — h) 



4F3(l,l,l,l;a+l,6,x;l) 
4F3(l,l,l,l;a,6+l,x;l) 
4^3(1, 1,1,1; a, 6,x + l;l) 



(30) 



with a, 6 G N\{0, 1} and a^h. 



B Corrections from virtual and real emission 



The missing piece for the virtual corrections to the hard part of the partonic cross section 
for the [8a] configuration (see Eq. (pT]) ) is 

3ai(r) (r2-5)6i(r) {r-?,){r + lfh2{r) (4r - 13)63 (r) 

•^qq^&ay) — — — — TTTT- H — — ^ 



4(r-l) 12(r-l)2 



24(r- 1)2 



12 



2(r2 - 2r + 5)64(r) lQrh^{r) ^ Sbeir) ^ 2(r - l)6^(r) 



3(r2-l) 3(r2-l) 3(r + 1) ' 3 

(r + l)2(r2 - 2r + 5)ci(r) (r^ - 5r2 + llr - 15)c2(r) 



24(r-l)2 12(r-l)2 
3(r + l)(r2 - 6r + 17)c3(r) 3(r + l)(r2 - 2r + 5)c^{r) 



8(r-l) 4(r-l) 
(3r2 - 4r - 17) ln(2) ISr^ + 14r2 + 23r - 101 



4(r - 1) 



12(r2 - 1) 



+nf 



(r — 5)ai(r) rbi{r) {r + iyb2{r) 
6(r-l) 3~ ^ 6(r- 1) 

+ (r - l)Oi(r) + 



18(r - 1) 



(31) 



40 



where the scalar one-, two- and three-point functions are 



Re{i?oK;m|,0)} 
Re {5o(4m?;m?,m?)} 
5o(0; m|,m2) 
5o(-m?;m?,0) 

Re{Bo{ml;mlO)} 
lRe{B',{ml;mlO)} 



m 



m 
m 



2 2 2\ 1 



Re{Co(4m2,0,0; 
Re {Co(0, -m?,m?;m?,m?,0)} 
m?Co(0, -m?, m?; m?, m?, 0) 

C4(r) 

? Re {Co(0, -m?, m?; m?, m?, 0)} 



m? A + m? ai (r) 
A + 6i (r 
A + 62 (r 
A + 63 (r 
A + 64 (r 
A + 65 (r 
A + fee (r 

&2 {r) , 
ci (r) , 
C2 (?^) , 

m?Co(4m?, 0, ; m?, m?, m?) 



1 



r + 1 



1 + ln 



r + 1 



C3 (?^) , 
- In (r 



(32) 



with the conventions of [69] and A = l/^uv — 7^ + [47r/i|./(2mg)^]. The remaining 
functions complete the real corrections 



(r - l)(r + ifz^ - 2(r + l)(19r + 35)^^ 
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9r(r + I'fz^ 
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z(r+l) 
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C Benchmark scenarios 



This Appendix contains detailed information about the scenarios (a)-(q), (X) and (Y). 
The values of the msugra parameters defining the scenarios are listed in Tab. [5l In Tabs. [6] 
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benchmark 


SPS 


mo [GeV] 


mi 12 [GeV] 


[GeV] 


tan(/3) 


sign(^) 


point 


scenario 












(a) 


SPSla-point 


100 


250 


-100 


10 


1 


(b) 


SPSla-slopcl 


80 


200 


-80 


10 


1 


(c) 


SPSla-slope2 


60 


150 


-60 


10 


1 


(d) 


SPSla-slopc3 


120 


300 


-120 


10 


1 


(e) 


SPSla-slope4 


140 


350 


-140 


10 


1 


(f) 


SPSla-slopc5 


160 


400 


-160 


10 


1 


(g) 


SPSla-slope6 


180 


450 


-180 


10 


1 


(h) 


SPSla-slope? 


200 


500 


-200 


10 


1 


(i) 


SPSlb-point 


200 


400 





30 


1 


(j) 


SPS2-point 


1450 


300 





10 


1 


(k) 


SPS2-slopel 


1250 


200 





10 


1 


(1) 


SPS2-slope2 


1050 


100 





10 


1 


(m) 


SPS2-slope3 


1650 


400 





10 


1 


(n) 


SPS2-slope4 


1850 


500 





10 


1 


(o) 


SPSS-point 


90 


400 





10 


1 


(P) 


SPS4-point 


400 


300 





50 


1 


(q) 


SPS5-point 


150 


300 


-1000 


5 


1 


(X) 




900 


550 





10 


1 


(Y) 




400 


600 





10 


1 



Table 5: Initial parameters for the SPS scenarios and scenarios (X) and (Y) to obtain the 
benchmark points. 



and [7] we list the values for the squark masses as provided by SuSpect [38]. Note that for 
our analysis the averaged values as provided in Tab. |3]are used. 
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benchmark 
point 


mj^ [GeV] 


m^^ [GeV] 


niu^ [GcV] 


mu^ [GcV] 


m,^ [GeV] 


nis^ [GeV] 


(a) 


567.77 


545.62 


562.26 


545.89 


567.77 


545.62 


(b) 


464.72 


446.20 


457.94 


445.62 


464.72 


446.20 


(c) 


360.24 


344.94 


351.37 


343.26 


360.24 


344.94 


(d) 


669.23 


643.23 


664.59 


644.23 


669.23 


643.23 


(e) 


769.54 


739.52 


765.52 


741.20 


769.54 


739.52 


(f) 


868.90 


834.78 


865.36 


837.09 


868.90 


834.78 


(g) 


967.33 


929.03 


964.15 


931.94 


967.33 


929.03 


(h) 


1065.22 


1022.71 


1062.35 


1026.20 


1065.23 


1022.71 


(i) 


876.84 


843.16 


873.27 


845.41 


876.84 


843.16 


(j) 


1556.70 


1552.46 


1554.76 


1552.07 


1556.70 


1552.46 


(k) 


1300.76 


1300.38 


1298.41 


1299.46 


1300.76 


1300.38 


(1) 


1055.59 


1058.08 


1052.65 


1056.68 


1055.59 


1058.08 


(m) 


1817.49 


1808.85 


1815.83 


1809.03 


1817.49 


1808.85 


(n) 


2080.64 


2067.30 


2079.20 


2068.06 


2080.64 


2067.30 


(o) 


858.66 


824.01 


855.07 


826.36 


858.66 


824.01 


(P) 


764.91 


743.68 


760.88 


744.40 


764.91 


743.68 


(q) 


678.15 


652.40 


673.86 


653.45 


678.15 


652.40 


(X) 


1437.13 


1403.17 


1435.06 


1406.01 


1437.13 


1403.17 


(Y) 


1296.51 


1247.59 


1294.17 


1252.01 


1296.51 


1247.59 



Table 6: d, u, s masses for the various scenarios. The numbers have been obtained with 
SuSpect [38] 
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benchmark 
point 


ms^ [GeV] 


ms^ [GeV] 


^i, [GeV] 


^l, [GeV] 


mi, [GeV] 


lui^ [GeV] 


(a) 


562.26 


545.89 


516.91 


546.24 


399.73 


586.53 


(b) 


457.94 


445.62 


421.76 


448.10 


319.81 


499.59 


(c) 


351.37 


343.26 


325.09 


348.45 


242.07 


413.19 


(d) 


664.59 


644.23 


610.48 


642.77 


479.61 


673.26 


(e) 


765.52 


741.20 


702.91 


738.14 


559.00 


759.85 


(f) 


865.36 


837.09 


794.42 


832.56 


637.74 


846.32 


(g) 


964.15 


931.94 


885.04 


926.03 


715.69 


932.52 


(h) 


1062.35 


1026.20 


975.15 


1018.98 


793.11 


1018.75 


(i) 


873.27 


845.41 


777.67 


826.04 


661.71 


839.81 


(j) 


1554.76 


1552.07 


1298.73 


1539.74 


970.87 


1307.45 


(k) 


1298.41 


1299.46 


1073.69 


1289.26 


789.69 


1083.38 


(1) 


1052.65 


1056.68 


860.37 


1048.57 


622.94 


870.89 


(m) 


1815.83 


1809.03 


1528.82 


1794.54 


1157.84 


1536.48 


(n) 


2079.20 


2068.06 


1761.18 


2051.40 


1347.10 


1767.83 


(o) 


855.07 


826.36 


792.13 


822.99 


649.29 


842.40 


(P) 


760.88 


744.40 


617.79 


685.28 


546.52 


696.19 


(q) 


673.86 


653.45 


561.85 


650.48 


248.72 


649.57 


(X) 


1435.06 


1406.01 


1279.28 


1394.96 


1035.19 


1300.47 


(Y) 


1294.17 


1252.01 


1189.74 


1242.61 


9878.50 


1221.42 



Table 7: c, 6, t masses for the various scenarios. The numbers have been obtained with 
SuSpect [38] 
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